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^ ! ABSTRACT 
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We report results from the observations of the well studied TeV blazar Mrk 421 
with the Swift and the Suzaku satellites in December 2008. During the observation, 
Mrk 421 was found in a relatively low activity state, with the corresponding 2 — lOkeV 
r~| ' flux of 3 x 10~ 10 ergs -1 cm -2 . For the purpose of robust constraining the UV-to-X-ray 

emission continuum we selected only the data corresponding to truly simultaneous time 
intervals between Swift and Suzaku, allowing us to obtain a good-quality, broad-band 

■ spectrum despite a modest length (0.6ksec) exposure. We analyzed the spectrum with 
the parametric forward-fitting SYNCHROTRON model implemented in XSPEC assuming 
two different representations of the underlying electron energy distribution, both well 
motivated by the current particle acceleration models: a power-law distribution above 

^sO ! the minimum energy 7 m i n with an exponential cutoff at the maximum energy 7 maX ; 

and a modified ultra-relativistic Maxwellian with an equilibrium energy 7 eq . We found 
that the latter implies unlikely physical conditions within the blazar zone of Mrk 421. 

0^ • On the other hand, the exponentially moderated power-law electron distribution gives 

two possible sets of the model parameters: (i) flat spectrum dN' e /d r y oc 7 -1 - 91 with 

7— l ■ low minimum electron energy 7 m i n < 10 3 , and (ii) steep spectrum oc 7~ 2 - 77 with high 

• • , . 

. minimum electron energy 7 m i n — 2 x 10 . We discuss different interpretations of both 

possibilities in the context of a diffusive acceleration of electrons at relativistic, sub- or 

superluminal shocks. We also comment on how exactly the 7-ray data can be used to 

discriminate between the proposed different scenarios. 



department of High Energy Astrophysics, Institute of Space and Astronautical Science (ISAS), Japan Aerospace 
Exploration Agency (JAXA), 3-1-1 Yoshinodai, Sagamihara, 229-8510 Japan 

2 Department of Physics, University of Tokyo, Hongo 7-3-1, Bunkyo, 113-0033 Japan 

3 Astronomical Observatory, Jagiellonian University, ul. Orla 171, Krakow 30-244, Poland 

4 Stanford Linear Accelerator Center/KIPAC, 2575 Sand Hill Road, Menlo Park, CA 94025, USA 

5 Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo, 169-8555 
Japan 



- 2 - 



Subject headings: acceleration of particles — radiation mechanisms: non-thermal - 
galaxies: active — BL Lacertae objects: individual (Mrk 421) — galaxies: jets - 
X-rays: galaxies 



1. Introduction 



Blazars constitute a sub-class of radio-loud, jet-dominated Active Galactic Nuclei (AGN) char- 
acterized by a broad-band non-thermal emission detected from radio to high energy 7-ray band, 
strong variability on time scales ranging from hours up to years, two distinct (low- and high-energy) 
spectral humps in the v-vF v representation, and also by a substantial polarization in the radio and 
optical photon energy ranges. The low-energy spectral component of blazar sources, peaking at 
IR-to-X-ray frequencies, is widely believed to be due to a Doppler-boosted synchrotron emission 
of ultra-relativistic electrons accelerated in the innerm ost parts of relativistic and oriented close 
to our line of sight jets containing magnetized plasma (|Urrv &: Padovanilll995l ). The high-energy 
spectral component of low-power blazars of the BL Lacertae type (hereafter BL Lacs) extends up 
to a few or even tens of TeV, and is most commonly considered as originating via inverse-Compton 
scattering by ultra-relativistic electrons of synchrotron photons produced by the same electrons 
("synchrotron self-Compton process," hereafter SSC). This 7-ray emission, which only recently can 
be studied in detail with the new generation of ground-based Imaging Atmospheric Cherenkov 
Telescopes (IACTs) and Fermi satellite, exhibits particularly dramatic v ariability, observed on the 
timescales as short as minutes (lAharonian et alJl2007l ; Albert et al.ll2007l ) . This indicates extremely 
efficient and rapid particle acceleration processes taking place in compact emission regions. Our 
understanding of such acceleration processes is, on the other hand, still insufficient, despite a sig- 
nificant theoretical progress made in this field during the last years. 



Markarian 421 (Mrk 421; z = 0.03), the first established extragalactic TeV source (jPunch et al 



1992), is one of the most famous and nearby "high frequency peaked" BL Lacs. Since its discovery, 
it has been routinely observed by several IACTs at very high energy 7-rays and, more importantly, it 
has been targeted by many intensive multi- wavelength campaigns (which, in the past, were limited 
in their temporal and spectral coverage, and notably, often were triggered by a particularly high 
activity state of the source). These resulted in the discovery of a significant correlation between 
the hard X-ray and TeV fluxes, thus supporting the SSC origin of its high-energy emission (e.g. 



Inoue fc Takaharal 1 19961 : iTakahashi et al. 



19961 : iFossati et al.l 120081 ) . In a framework of such an 



interpretation, the energy distribution of the radiating electrons is required to be of a broad, non- 
thermal form, and the power-law spectral shape is the most commonly invoked approximation for it. 
The main justification for this approximation is that the diffusive shock acceleration, the "lst-order 
Fermi process" — known to generate power-law spec tra of ultra-relativistic par t icles — is likely 



to play a major role in accelerating jet electrons (e.g. iMastichiadis fc Kirld 119971 : ITakahashi et al 



- 3 - 



200Q : ISpada et al.ll200ll ). Still, as pointed out by several authors, alternative forms of the electron 
energy distribution in blazar jets, namely "ultra-relativistic Maxwellian" expected to be produced 
in the case of the efficient turb ulent ("2nd-order Fermi") acceleration, may als o account for bulk of 
the observational findings (e.g. ISauge fc Henril 12004 ; iKatarzvhski et al.l l2006bl : iGiebels et al J 120071 : 
Ushio et al. I bood ). In order to determine which acceleration process is likely to be dominant, it 
is necessary to reconstruct the underlying electron spectrum from the observed emission spectra. 
The X-ray observations are crucial in this respect. 



The Suzaku X-ray observatory (jMitsuda et al.l 120071 ) has the capability of high quality obser- 
vations of astrophysical obj ects from 0.5 up to 3 00 keV, owing to the excellent sensitivity of its 
Hard X-ray Detector (HXD: iTakahashi et al.ll2007l ). As such, it is the most suit able instrument fo r 
analyzing the synchrotron continua of high-energy peaked BL Lacs. Recently, lUshio et al.l (|2009) 
demonstrated the excellent quality of the Suzaku data by constructing high-photon-statistics time- 
resolved 0.6 — 60keV spectra of Mrk 421 for exposure times as short as a few ksec. These were 
analyzed by the parametric forward-fitting SYNCHROTRON model (rather than being fitted by some 
assumed spectral shape of the synchrotron continuum), to investigate more directly the electron 
energy distribution. Here, we proceed further in this direction, by fitting the broad-band spectrum 
of Mrk 421 with different forms of the electron distribution discussed in the literature. In addition 
to the new Suzaku data (collected d uring December 200 8), however, we also utilize the simultaneous 
observations by the Swift satellite (IGehrels et al.ll2004l ). and in particular the Ultra- Violet/Optical 
Telescope (UVOT: iRoming et al.ll2005l ) onboard Swift. §[2] reports in more detail on the analyzed 
datasets. In §(3j we present the results of the applied spectral analysis. The implications of our 
results towards the understanding of blazar sources are discussed and summarized in §[4] and §0 
Throughout this analysis, we use the XSPEC version 12.5.0ac and the CALDB of version 2009.04.03. 



2. Observations and Data Reduction 



2.1. Suzaku 



Suzaku observations were performed with the X-ray Imaging Spectrometer fXI S:lKovama et al 



2007 



2007) in the 0.3 — 12keV band, and also with the Hard X-ray Detector (HXD; ITakahashi et al 



Kokubun et al.ll2007l) in the 10 — 300 keV band. The XIS cameras, located at the foci of the 



X-ray telescopes (XRTs; iSerlemitsos et al.l 120071 ) . include one back-illuminated CCD camera (BI) 
and two front-illuminated CCD cameras (FIs). The HXD is a non-imaging detector system which 
consists of silicon PIN diodes capable of observing in the 10 — 70 keV band and the GSO crystal 
scintillators in the 40 — 300 keV band. In the following analysis, we do not use the GSO data, 
because Mrk 421 was not sufficiently bright above 40keV during the performed observations to be 
significantly detected by the GSO. 

We observed Mrk 421 with Suzaku for almost three days from MJD 54803.76 to MJD 54805.89 
(seqnum=703020010). All of the XISs are operated with 1/4 window option in order to reduce 
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possible pile- up effects. Both XIS and HXD data are reprocessed from "uncleaned event files" 
with the CALDB of 2009-04-07 which corresponds to rev2.2 of the official pipeline process. The 
"unscreened event files" of XISs are reduced with the standard selection criteria requiring, in 
addition, the rigidity larger than 6GV/c. The XIS events are next extracted from a circular region 
with a radius of 3' centered on the image peak. The background region is extracted from the source 
free region having a similar area. Moreover, we extract the image center of BI sensor with a radius 
of 0.2' because it weakly suffers from pile-up judging from its count rate. In our spectral fitting, 
we derive the response matrices (RMF) and the effective area file s for the XRTs (AR F) for each 
XIS sensor by using xisrmfgen and xissimarf gen, respectively (jlshisaki et al.l 120071 ). Since the 
two FI CCDs have almost identical performance, we co-add their data and the corresponding RMF 
and ARF. For the HXD data, we exclude events during South Atlantic Anomaly (SAA) passages, 
Earth occultation and those with cutoff rigidity less than 6GV/c. 



2.2. Swift 



The three-day-long Suzaku observation of Mrk 421 is well complemented by two Swift observa- 
tions, involving nine short exposures (~ 200s each; seqnum=00030352102 and 00030352103). The 
Swift satellite carr ies three sets of instruments: the Burst Alert Telescope (BAT ; 15 — 150 keV; 



Burrows et al.ll2005l ). and the Ultra 



200 



5]). Hereafter, we only an- 



Gehrels et al.ll2004h . the X-ray Telescope (XRT; 0. 3-10keV; 
Violet/Optical Telescope (UVOT; 170 - 650 nm; iRoming et al. 
alyze the XRT and the UVOT data. All the data are taken from the NASA/ Swift database 
(http://heasarc.gsfc.nasa.gov/FTP/swift/) and are processed at the Swift Data Center. 

The Swift XRT data were all taken in the Windowed Timing mode, with the central 8 arcmin 
field of view covered, and one dimensional imaging preserved. Therefore, we extract the rectangular 
region of 0.50 x 3.0 arcsec 2 centered on the peak position of the source. The background is selected 
in the source free region with rectangular shape. Photons with grade — 2 are selected and 
rebinned so that each bin has sufficient photon statistics. The auxiliary response file is created by 
the task xrtmkarf and the standard response file swxwt0to2s6_20070901v011 . rmf . For UVOT 
observations, we utilize the three ultraviolet filters in the imaging mode, with the effective central 
wavelength of 260.0 nm, 224.6 nm and 192.8 nm for UVW1, UVM2 and UVW2 filters, respectively. 
The source aperture sizes are all chosen to correspond to those used to determined the UVOT 
zero points, namely 5" for the ultraviolet filters. A 15" background region is extracted from the 
"ghost wing" of nearby bright source 51 UMa. A comparison of results using seve ral different 



backg round regions suggests that the corresponding uncertainty is at most 0.01 mag (jPoole et al 



2008), which is smaller than the systematic uncertainty indicated by the UVOT team. Next, we 
calculate the photometric spectra so that they can be fitted in XSPEC by the task uvot2pha (vl.5) 
using the response matrices swu[wl ,m2 , w2] J20041120vl04 . rsp. Note that the columns "E_MIN" 
and "E_MAX" of these response are manually set to the FWHM of the filter in order to improve 
their appearance. 
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3. Analysis 

3.1. The SYNCHROTRON Model 

The broad bandpass and high-sensitivity of the Suzaku satellite enables spectral fits for the 
observed X-ray emission of celestial sources with synchrotron spectra originating from a given 
(assumed) shape of the electron energy distr ibution (EED), rather than by the as sumed spectral 
form of the non-thermal emission continuum (jTanaka et al.ll2008l ; lUshio et al.ll2009l ). We note that 
in the past the analogous "parametric forward-fitting" technique has been applied in analyzing 
X-ray emission of solar flares, believed to be produce d by bremsstrahlung of non-thermal elect rons 
within the magnetic loops above the solar surface (see Park et al. 1997 ; Petrosian k, ChenlholO , and 
re ferences ther e in) . In this paper, we continue the studies of the X-ray emission of Mrk 421 initiated 



m lUshio et al.l (|2009i ). including the ultraviolet observations by UVOT onboard Swift satellite in 
addition to the analyzed Suzaku dataset. Below we briefly describe the SYNCHROTRON model serving 
this function, as implemented in XSPEC. 

In our SYNCHROTRON model, we assume two different forms for the EED, both well motivated 
by the particle acceleration models discussed in the literature (see next section), namely a power- 
law (PL) form above the minimum energy 7 m ; n with an exponential cutoff at the maximum energy 
Tmax; and a modified ultra-relativistic Maxwellian (UM) with the equilibrium energy 7 eq : 



d'y 



exp 



7max 



NqJ 2 exp 



7 

7eq 



for 



for 



1 <7n 

1 < 7 



< 



7 



(PL) 
(UM) 



(1) 



where 7 = E' e /m e c 2 is the electron Lorentz factor, Nq is the normalization of the EED, s is 
the electron energy index, and a represents the modulation of the exponential cut-off. Here and 
below the primes denote the quantities measured in the emitting plasma rest frame, except of 
electron Lorentz factor 7, which is always measured in the this frame. The considered model PL 
is thus characterized by three free parameters (s, 7 m i n , and 7 max ), while the model UM by two 
free parameters only (a and 7 eq ; since with the extremely flat low-energy spectrum oc 7 2 the exact 
location of 7 m i n does not affect the modeled synchrotron continuum). For these, we consider the 
following "expected" ranges: 1.5 < s < 3.5, 1 < 7 m i n < 7 max , 10 4 ' 5 < 7max < 10 6 ' 5 , 0.1 < a < 2, and 
7eq < 10 5,0 . We note that this is the first time when the two aforementioned representations 
of the EED are tested against a single, broad-band dataset. 

Having specified the form of the EED with the assumed isotropic pitch-angle distribution, the 
jet comoving synchrotron emissivity j' syn (e') is described by: 

dN: F(i ]d -„ (2) 



where the function F(x) is defined as 

F(x) 



(3) 
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and where ^5/3 (£) is the modified Bessel function of the 5/3 order, the characteristic synchrotron 
photon energy e c is 

, 3heB' 2 , . 

4 = - A 1 , (4) 

and B' is the jet magnetic field intensity. In the case of a relativistic bulk motion of the emitting 
plasma, the observed energy of the synchrotron photon is e = S e', and the spectral emissivity 
transforms as j sy n(e) = <5 2 jsyn( e ')> where 5 is the corresponding Doppler facto4j. In the following 
analysis we adopt the product of B'5 = 2 (where B' is expressed in units of Gaus s) for the blazar 



emissi on zone in Mrk 421, which is the value claimed previously for this source bv lTakahashi et al 



(|2000l ) based on detailed variability and spectral studies, when the object showed a similar flux to 
that measured by us. The validity of adopting B'5 = 2 (and an investigation of a range for values 
adopted for this parameter) is discussed later. Hence, the observed synchrotron flux at a given 
observed photon energy in our model depends solely on the parameters of the EED. 



3.2. Spectral Fitting 

Our observations took place during a relatively low activity state of Mrk 421, characterized 
by only a mild variability: the observed 2 — lOkeV flux varied between ~ 2.7 x 10~ 10 and ~ 
4.1 x 10~ 10 ergs^ 1 cm ~ 2 during the 3-day long Suzaku observations. In order to construct a 
truly simultaneous broad-band spectrum ranging from UV to hard X-ray frequencies, we have 
extracted the available data taken only during the overlapping time intervals for XRT, XIS and 
HXD/PIN instruments; this resulted in 0.6ksec of the final net exposure. Nevertheless, even 
with such a short exposure, we were still able to obtain a very good quality broad-band X-ray 
spectrum of Mrk 421 extending from 0.5 up to 50keV. The averaged X-ray spectrum constructed 
by accumulating all of the detected photons instead would not be significantly different than this 
truly simultaneous spectrum because, as noted above, Mrk 421 was relatively steady during the 
observations, and no significant spectral changes have been noted in the total accumulated Suzaku 
dataset. The conventional fit to the X-ray spectrum with a single power-law function modified 
solely by the Galactic absorption describes the data poorly (xt = 3.47 with 56 dof), and indicates 
a significant curvature of the X-ray spectrum. A broken power-law function provides a better 
representation of the data (xt = 1-06 with 54 dof), returning the low- and high-energy photon 
indices of Ti ow = 2.33 ± 0.01 and Thigh = 2.71 ± 0.04, respectively, with the break photon energy 
of ebrk = 2.82 ± 0.17keV. The corresponding 2 — lOkeV ener gy flux of Mrk 421 is th en S2-iokeV — 



2.83 x 10 ergs cm , again, similar to that reported bv lTakahashi et al.l (2000) 



The constructed dataset is then analyzed using the SYNCHROTRON model. The top panel of 
Figure Q] presents the UV-to-X-ray spectral energy distribution of Mrk 421 fitted with three different 



1 These transformations, together with the transformation of the emitting volume V = 5 V' as appropriate for a 
homogeneous moving source assumed here, give the energy flux transformation [eSe] = 5 4 [e S' e i], since [eS c ] oc [ej e ] xV. 
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representations of the EED as described below, while the bottom panels show the residuals between 
the data and the models. The spectral fitting is performed jointly in the UVOT band and the 
0.5—50 keV band; here, we ignore the data within the 1.6—2. keV and the 10—15 keV photon energy 



range s, because of the systematic calibration uncertainties (jSerlemitsos et al.ll2007l ; iFukazawa et al 
The extinction and the absorption of the UV/X-ray photons in our Galaxy is taken into 



account by means of multiplicative constants based on the analysis by iPeil (119921 see Table [JJ . 
Such a procedure w as applied previously in the analysis of the spectra of Gamma-Ray Bursts (e.g. 



Schadv et al.l 120071) . An alternative analysis of the UV data following the procedure presented in 



Kataoka et al.l (|2008l ) returns the fluxes consistent within a few percent with the ones obtained by 



our method. 

In the first representation of the fitted EED — denoted below as the case "a" — we consider 
the power-law (PL) model with fixed minimum electron Lorentz factor 7 m i n = 1. The results of 
the model fitting are summarized in Table [T] and shown in Figure [TJ The emerging parameters 



are s = 1.91 t^ra and log(7 max ) = 5.09I^t, with the quality of the fit xt = 1-86 (for 54dof). 
The inset in Figure [JJ shows that these two parameters (s and 7max) are weakly correlated with 
each other in a sense that the harder the electron distribution gets, the lower maximum energy of 
electrons is required to fit the data. However, the remaining residuals around 0.5 keV indicate that 
the case "a" may not be the best representation of the data. We emphasize again that here we 
probe directly the parameters of the EED, unlike the previous studies in which the observational 
data were fitted by some assumed parametric form o f the emission continuu m (either a power-law, 
a broken power-law, or a log-parabolic form; see e.g.. iTramacere et al.ll2009i ). In addition, the EED 
parameters are determined robustly, because the UV data are taken into account in addition to the 
X-ray data. Without truly simultaneous UV data, the model fits would have a larger degeneracy. 



In the second representation of the fitted EED (hereafter case "b") we set free all the three 
parameters of the power-law model PL. The resulting values are then s — } ~~~ <li " 



4.32+0^ and log(7 max ) = 5.36+^ with Xv = 1-12 (for 53dof). The quality of the fit is now 
significantly improved, as shown in Figure [TJ Finally, in the third representation of EED denoted 
as the case "c", we consider the modified ultra-relativistic Maxwellian form (model UM). In this 
case, which represents the data well (xt = 1-20 with 53dof; see again Tabled] and Figure [JJ, the 



model parameters are log(7 eq ) = 1.41_qq3 and 1/a 



5 14+0.06 

■J- 1 -0.09- 



3.3. Related Uncertainties 

The model parameters emerging from the analysis described above may depend to some level 
on the assumed value of the Galactic column density in the direction of Mrk 421, which affects 
the low energy X-ray data. Therefore, it is necessary to investigate how robust are the presented 
model fits by varying Nn within the allowed range and repeating the analysis. In the case of the PL 
models (cases "a" and "b"), Figure [2] shows the resulting dependence of the electron energy index 
s on the level of the Galactic absorption, with the extinction parameter E(B-V) varied accordingly 
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Table 1: The fitting results with the SYNCHROTRON model 3 



model 


case 




model parameters 






2 - lOkeV flux 


xl (dof) 






log (7 m i n ) 


l0g(7max/eq) 


s 


or 1/a 


52-iokcV [crgcm~ 2 s -1 ] 




PL 


"a" 


0.0 (fixed) 


5.09 ±0.01 


s = 


1 qi+0.01 


2.61 x 10~ 10 


1.86 (54) 


PL 


"b" 


4.32 ±0.03 


c oc+0.06 
°- OD -0.05 


s = 


9 77 +0.13 
z - ' ' -0.12 


2.65 x 10~ 10 


1.12 (53) 


UM 


"c" 


0.0 (fixed) 


1 41+0.06 


1/a = 


- 5 1 4+O.O6 


2.65 x 10~ 10 


1.20 (54) 



"Obtained with the multiplic ative constants co nst=UVOT:XRT:XIS:PIN =1.00:1.07:1.00:1.13, E(B -V)=0.015, 
R V =3.1, and z=0.0 for zdust (jSchlegel et al.lll998h , and N H =0.0161 for wabs (|Lockman fc Savage|[l99i ). 



(assuming linear scaling with Ah). In the fig ure, the red box indicates the allowed range for Ah 



estimated from the historical HI observations (IDickey &: Lockmanlll990l : lLockman Savage 



1995 



Kalberla et al.ll2005l ). As is shown there, in the case "a" (filled circles in the plot) the slope of the 
EED is almost independent on the particular value of the absorption/extinction coefficients, but 
the acceptance of the fit gets significantly worse as Ah becomes larger. This behavior is easy to 
understand, since with the fixed minimum electron energy 7 m ; n = 1 the remaining two parameters of 
the power-law model are determined by the broad-band UV-to-X-ray emission continuum, almost 
irrespectively on the low-energy segment of the X-ray dataset (< 1.0 keV) which might be most 
affected by increased Ah and which, as such, determines only the quality of the fit. The same 
is true for the UM model (case "c", not included in Figure [2]); namely, changes in the hydrogen 
column density (and, correspondingly, the extinction level) do not change significantly the emerging 
parameters of this model, but only affect the quality of the fit. We note that with the assumed 
linear scaling of E(B-V) with Ah the fluxes in the UVOT band vary much less substantially (i.e., 
within 10%) than the soft X-ray flux. 

On the other hand, in the case of the PL model with all three parameters set free (case "b", 
open circles in Figure [2]), the emerging electron energy index correlates with Ah, while the quality 
of the fit remains very good. As a result, within the allowed range of the Galactic absorption the 
electron index in the case "c" reads roughly as s ~ 2.5 — 3.0 (compare it with Tabled]). Again, the 
observed dependence is not difficult to understand. This is because in the case "b" the minimum 
electron Lorentz factor is so large that the observed synchrotron emission in the UV frequency 
range is produced almost entirely by the lowest-energy electrons. In particular, with the assumed 
B'5 = 2 the observed characteristic frequency of the synchrotron photons emitted by the electrons 
with 7 m i n = 2 x 10 4 is e c ~ 15 eV. In other words, the slope of the electron energy energy distribution 
in the case "b" is determined not predominantly by the broad-band UV-to-X-ray emission, as in 
the cases "a" and "c", but solely by the X-ray continuum. As such, the parameter s in the case 
"b" is sensitive to the low-energy X-ray photons, and therefore to the particular value of Ah- 

Another concern arises from the fact that in our analysis the product of the magnetic field 
intensity and the jet Doppler factor has been fixed at the particular value B'5 = 2. There have been 
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numerous efforts to determine uniquely the physical parameters of blazar emissio n region in TeV- 



detected BL Lacs, in cluding B' and 5 , in ad dition to the aforementioned work bv lTakahashi et al 



(|2000l ). For example. iTavecchio et al.l (| 19981 ). who attempted to utilize the multi- wavelength data 
self-c onsistently, claimed 5 = 25 and B' = 0.2 G for the flari n g stat e of Mrk 421 (see lTakahashi et al 
19961 ). This would lead to B'6 = 5. Similarly, lAlbert et al.l feooj ) obtained B' = 0.2 G and 5 = 15 



resulting in B'5 = 3. Although those values are slightly different from the one assumed in the our 
paper, they are of the same order of magnitude. In order to investigate in detail how the related 
uncertainty regarding the product B'5 affects the presented results, we repeated the analysis for 
B'8 = 1,5, and 10. We found that within the range B'5 = 1 — 10, the resulting curvature of the 
EED (parameters s or a) does not significantly change when compared to B'5 = 2, but only that 
the critical electron energies (7 m i n , 7max and 7 eq ) scale as oc \ j\J B' 5, as expected (see equation H]). 
Therefore, we can safely conclude that the uncertainty regarding B'5 does not affect significantly 
the model results. 



4. Discussion: Insights into the Energy Distribution of Radiating Particles 

In the previous sections we showed that the UV-to-X-ray synchrotron continuum of Mrk 421 
observed in December 2008 may be fitted almost equally well by three different representations of 
the EED, each well motivated on theoretical grounds, yet each implying different physical conditions 
and different particle acceleration processes operating within the blazar emission zone of Mrk 421. 
Here we discuss in this context and in more detail the results obtained by us. In the following, we 
assume that the inferred shapes of the EED may be identified with the steady-state spectra formed 
within a homogeneous acceleration/emission region. It is not clear (or it is even questionable), 
however, whether this particular assumption — even though being widely applied in blazar modeling 
- is a good approximation of a realistic situation, and we comment on this issue in § 14.41 below. 



4.1. Ultra-relativistic Maxwellian 



We begin with the case "c" analyzed above, namely with the ultra-relativistic Maxwellian rep- 
resentation of the energy distribution of the radiating electrons. In general, it was shown previously 
by various authors that the Maxwellian spectra oc r y 2 exp[— 7/7eq] are expected to form in the case 
of an efficient continuous acceleration of relativistic particles undergoing radiative (synchrotron- 
type) losses obeying r^ ad oc 7 -1 and a very inefficient escape from the acceleration/emission region 



(Schlickeiser 



1984 



Aharonian et al.lll986l : iPark fe Petrosianlll995l ). In particular, turbulent ("2nd- 
order", or "stochastic") acceleration of the Fermi-type — characterized by the acceleration timescale 
independent on the particle energy, r^ cc oc const — was the most commonly considered example of 
the continuous energization of ultra-relativistic particles, with t he equilibrium energy 7 P q d efined 



via the condition r' acc = T r ' ad (7 e q). This has been generalized by lStawarz fc Petrosianl (|2008l ). who 



showed that for the power-law scalings t^ cc oc 7 p and r' d oc j q , the steady-state solution to the 
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appropriate equation governing energy evolution of the radiating particles (assuming, again, very 
inefficient particle escape from the system) is of the modified ultra-relativistic Maxwellian form 
oc 7 2 exp[— (1/a) (7/7eq) a ], where a = p — q. 

Our analysis indicates that in order to fit the UV-to-X-ray spectrum of Mrk 421 with the 
modified Maxwellian form of the EED, the equilibrium energy has to be very small, 7 cq ~ 30, 
and the spectral curvature has to be very smooth, a ~ 0.2. Even though it is formally possible 
to relate these requirement to a set of particular physical conditions within the emission zone 
of Mrk 421, the derived parameters appear unrealistic. For example, anticipating the dominant 
inverse-Compton cooling to take place in the Klein-Nishina (KN) regime, i.e. q ~ 0.5, one would 
have to assume p ~ 0.7 in order to obtain a ~ 0.2. This requirement, interestingly, would then be 
consistent with the stochastic acceleration of the radiating particles by the turbulence characterized 
by the energy spectrum similar to the Kreichnan on^E However, the anticipation that relatively 
low-energy electrons within the blazar emission zone of Mrk 421 cool predominantly via the inverse- 
Compton emission in the KN regime rather than via synchrotron emission is highly questionable. 
For the case where the synchrotron emission dominates (q = — 1), on the other hand, no commonly 
acceptable value of p can be found to satisfy a ~ 0.2. Hence, we consider the case of the UV- 
to-X-ray spectrum of Mrk 421 fitted by the modified Maxwellian form of EED as only a formal 
possibility only, involving somehow unphysical (or, at least, "non-standard") conditions. 



4.2. Shock-originated Power-law 



Next we discuss the power-law models for the EED, focusing first on the case "b". In this case, 
the energy spectrum of radiating electrons is characterized by a steep slope s — 2.77 within a very 
(uncomfortably?) narrow energy range 2 x 10 4 < 7 < 2 x 10 5 (preceded by a sharp low-energy cut-off 
or eventually a very hard low-energy tail). At the first glance, such a representation of the EED may 
seem as artificial as the on e discussed in the p r evious section, but we note here that it is in fact very 
similar to that invoked by iKatarzvhski et al.l (|2006al ) for the case of a distant (z = 0.186) BL Lac 
object 1ES 1101-232 surprisingly detected at TeV energies, where the absorption of its TeV emission 
due to pair production on the extragalactic diffu se background light w as expected to provide large 
opacity at TeV energies (see also in this context iTavecchio et al.ll2009l . for the case of z = 0.14 BL 
Lac 1ES 0229+200). In addition, the high minimum electron energy and the steep electron spectral 
index (significantly different from the "canonical shock spectrum" s = 2.0) might in principle be 
both reconciled with the diffusive ("lst-order-Fermi") acceleration of the radiating electrons at a 



2 For the energy spectrum of the turbulence W(fc) oc k~ a , where k is the wavevector of the turbulent modes, 
stochastic accel eration timescale reads as r^ cc oc 7 2_<T . Hence, for the Kreichnan spectrum a = 3/2 one obtains 
(see Stawarz fc Petrosianl 20081 . and references therein). 



oc 7 



1/2 
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relativistic, perpendiculaj§, and proton-mediated shock. That is because under such conditions the 
inertia of protons carrying bulk of the jet energy (and thus shaping the shock structure) define 
the critical electron energy 7 cr ~ r s h m p /m e ~ 2 x 10 4 (r s h/10) (where T s h is the bulk Lorentz 
factor of the shock i n the upstream rest frame), below which no lst-order-Fermi acceleration is 
in principle allowed (jBegelman fc Kirk! Il990l ; iHoshino et al.l Il992l ) . It is thus interesting that the 
minimum electron energy obtained in the case "b" is in agreement with the expected value of 7 cr , 
but only in the case of ultrarelativistic shock velocity, T s h — 10. In addition, the steep power- 
law slope of the electrons with energies higher t han 7 nr is in agreement with those claimed for 
the d iffusive acceleration at superluminal shocks (jNiemiec Ostrowskill2004l ; ISironi &: Spitkovskv 
20091 ). assuring self-consistency of the interpretation. 



Even though similar arguments as above (based on the spectral forms of the EED derived 
from the broad-band modeling) were given p reviously in favor of a dominant role of perpendicular, 



proton-mediated shocks in AGN jets (e.g., IStawarz et al.l 120071 : ISikora et al.l |2009| ). the concern 



regarding the presented interpretation of the dataset analyzed in this paper is that large 7 m i n 
obtained in the case "b" may be an artifact of the analysis forcing the broad-band UV-to-X-ray 
data to be fitted by an EED of single power-law form. Indeed, as already noted in § 13.31 with the 
assumed B'S = 2 the observed characteristic frequency of the synchrotron photons emitted by the 
electrons with 7 m i n = 2 x 10 4 is e c ~ 15 eV, i.e. just above the UV frequency range. Thus, the 
obtained value for the minimum electron Lorentz factor in the case "b" may not reveal the nature 
of particle acceleration within the blazar emission zone of Mrk 421, but only the fact that the 
analyzed UV and X-ray data cannot be fitted together by a single power-law form of the EED. 



Besides the aforementioned caveat, the particular form of the EED discussed here would require 
a lot of fine tuning to ensure 1 <C 7 m in ^5 7max within the emission region, or in other words to 
prevent formation of a low-energy tail in the electron distribution (energies 7 < 7 m i n ) due to the 
radiative cooling of ultrarelativistic particles injected from the shock to the emission site. Note in 
this context that with a postulated efficient re-acceleration processes taking place downstream of 
the shock, which could potentially remove the cooling problem, one should not expect the radiating 
electrons to maintain their injection spectrum, but rather to form a continuum resembling modified 
Maxwellian discussed in § 14.11 Yet the possibility for preserving the "shock-like" EED with high 
7min far from the injection site cannot be rejected robustly. Clearly, such a controversy could in 
principle be resolved with the simultaneous Fermi observations, because the position of 7 m i n is the 
major factor determining the slope of the 7-ray continuum within the MeV-GeV photon energy 
range. In other words, Fermi/LAT data, if simultaneous with the UV-to-X-ray data analyzed with 
the SYNCHROTRON model, could help to discriminate between the case "b" and other analyzed cases, 
i.e. between the high or low minimum electron Lorentz factor. 



3 "Perpendicular/oblique" shocks are those with magnetic field lines oriented at right /large angles to the shock 
normal. Furthermore, if the projection of the upstream fluid velocity on the upstream magnetic field direction is 
smaller/larger than the light speed, the respective shocks are called "subluminal/super luminal" . 
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4.3. Flat Power-law as a Cooled Maxwellian? 



Finally, we consider the power-law model, case "a", in which the minimum electron Lorentz 
factor is fixed at 7 m ; n = 1, and the emerging power-law slope of the EED is s ~ 1.93. We note 
that because only the UV-to-X-ray data are fitted in our analysis, any other fixed value for 7min 
within the range 1 — 10 3 would not affect the determination of the parameter s, since the observed 
characteristic frequency of the electrons with Lorentz factors 7 < 10 3 is expected to be located 
below the UVOT frequency range (with the assumed B'5 = 2). Thus the physical case to discuss 



is a flat (s ~ 1.91) electron power-law between the energies 7 m i n < 10 3 and 



7n 



10 5 . There 



may be two interpretation for such. One is the diffusive shock acceleration at a (mildly) relativistic 
subluminal shock (without necessari ly significant proton content); such a configuration results in 
a ver y flat particle spectra (s < 2.0; iKirk Heavens! 1 19891 ; lOstrowskil 1 199 ll ; iBednarz &: Ostrowski 
19961 ). In a framework of this interpretation, the entire observed synchrotron continuum of Mrk 421 
would then be produced in a slow-cooling regime. 

The other possible physical interpretation for the case "a" discussed in this paper may be a 
radiatively re-processed (cooled) Maxwellian electron distribution. As is well known, with a steep 
(s > 2.0) power-law EED continuously injected into the emission region, synchrotron-type radiative 
cooling (r r ' ad oc 7 -1 ) results in further steepening of the injected spectrum (by As = 1.0) above 
some critical (break) energy. On the other hand, in the case of a very flat, high-energy EED 
injected to the emission region, radiative cooling leads to formation of a low-energy power-law tail 
oc 7~ 2 instead. For example, a radiatively re-processed ultra relativistic Maxwellian 7 2 exp[— 7/7^ 



would give a familiar 7 2 exp[— 7/7 eq ] form of the EED. Moreover, as pointed out bv lGiebels et al 



(|2007l ). additional non-negligible inverse-Compton/KN losses may significantly flatten such a power- 
law tail, leading to a formation of cooled power-law spectra with s < 2.0. We believe that this 
mechan ism is a viable optio n for the interpretation of the synchrotron emission of Mrk 421. And 
indeed, iGiebels et al.1 1)20071 ) applied the discussed scenario in the particular context of Mrk 421, 
showing that a single power-law with s = 1.7 — formed as a cooled Maxwellian — may account 
for the bro ad-band emiss i on of this blazar observed in March 2001. The value of the electron slope 
claimed bv lGiebels et al.l (|2007l ) is lower than the one found in our analysis but, interestingly, the 
average 2 - lOkeV flux of Mrk 421 in March 2001 (£2-10 keV - 13 - 4 x 10~ 10 ergcm~ 2 s -1 ), was 
about 5 times higher than during the exposure analyzed here. Such a behavior may be reconciled 
with the interpretation above if one assumes increasingly pronounced KN effects in cooling of the 
radiating electrons as the synchrotron flux increases. 

The above interpretation requires however a strong cooling regime for all the electrons down 
to energies of 7 ~ 10 3 . This, in turn, would imply relatively large emission region size R: by 
simply equating the dynamical timescale r dyn ~ R/c with the cooling (predominantly synchrotron 



by assumption) timescale r r ' ad 



— syn — 



67rm e c/o"T7-B , one may find very roughly 



R ~ 0.03 



VlO 3 / \0.5G 



pc. 



(5) 
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Since the corresponding variability timescale r var ~ R/c8 would then be of the order of days-weeks 
(for 7 = 10 3 , B' < 1 G, and the jet Doppler factor 5 < 10), any intraday flux variations could be 
hardly accounted by the discussed model. 

The interesting scenario outlined above can be easily and directly tested with the simultaneous 
UV, X-ray, and 7-ray data (obtained with the Fermi satellite and modern IACTs). That is because 
such a broad-band dataset would allow one to determine the ratio of the synchrotron and inverse- 
Compton cooling rates (or, more precisely, the relative contribution of the inverse-Compton/KN 
cooling), together with the parameters of the underlying EED using the SYNCHROTRON model. 
Note in this context that in the first interpretation of the case "a" — the one involving shock 
acceleration — no obvious correlation between the cooling rates ratio and the electron spectral 
index is expected, while in the second interpretation — the one involving radiatively re-processed 
Maxwellian component — a clear correlation between these two quantities is implied. We want to 
point out that there exists Fermi and MAGIC data during the 3-day long Suzaku observations. The 
analysis of these data set, as well as the physical interpretation of the overall SED in the context 
of the above presented model will be reported on a forthcoming publication. 



4.4. Steady Component? 



In the discussion above we have assumed that all the possible representations of the EED 
emerging from the analysis correspond to the steady-state spectra formed within a homogeneous 
acceleration/emission region of Mrk 421 jet. As already emphasized, this assumption may not be 
correct. And indeed, if the blazar emission zone is highly inhomogeneous, then the inferred electron 
spectra could instead be due to a superposition of multiple distributions formed within different 
sub-volumes of a system. In particular, the Maxwellian form discussed above — with the "prob- 
lematic" parameter a ~ 0.2 — may reflect some complex superposition of different Maxwellians 
with "physically acceptable" curvatures a > 0.5 or, alternatively, some multiply-broken power-law 
electron con tinuum. As such, i t may well correspond to the steady emission component claimed for 
Mrk 421 by lUshio et al.l (|2009l . see section 5 therein). Interestingly, the X-ray flux of this compo- 
nent and its spectral curvature (which we determined for the pu rpose of this pape r by fitting the 
UM form with the SYNCHROTRON model to the dataset analyzed in lUshio et all 12009 ) are both quite 
similar to those discussed here. 



Conclusions 



In this paper we have investigated which form of the electron energy distribution describes 
the best the UV-to-X-ray spectrum of Mrk 421 in its relatively low state (December 2008; the 
observed 2 — lOkeV flux ^-lokcV — 2.8 x 10~ 10 erg cm' 2 s ~ 1 )- I n our studie s we have applied the 
forward-fitting SYNCHROTRON model implemented in XSPEC (jUshio et al.ll2009 ); that is, the analyzed 
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dataset has been fitted directly with a synchrotron emission originating from a given (assumed) 
shape of the electron distribution, rather than by an assumed parametric form of the non-thermal 
emission continuum as typically presented in the literature. We have considered two different 
general representations of the electron spectrum, both motivated by the particle acceleration models 
discussed in the literature, namely a power-law form within some finite energy range modulated 
by an exponential cutoff at the maximum energy, and a modified ultra-relativistic Maxwellian. We 
have found that the latter one, which would correspond to the continuous (turbulent) acceleration 
of the radiating particles within the emission region with no particle escape, even though formally 
allowed (i.e. giving a reasonable fit to the data), implies rather artificial, and thus unlikely physical 
conditions within the blazar zone of Mrk 421. On the other hand, a power-law representation 
of the electron energy spectrum fits the data equally well. In this case, two possible sets of the 
model parameters have been found: (i) a flat-spectrum one oc 7~ L91 with low electron minimum 
energy 7 m i n < 10 3 , and (ii) a steep-spectrum one oc 7~ 2 - 77 with high electron minimum energy 
7min — 2 x 10 4 . We have discussed different interpretations for both, involving diffusive acceleration 
of the radiating electrons at relativistic, sub- or superluminal shocks, as well as radiatively re- 
processed electron distribution of the initial (injected) Maxwellian form. 

The presented technique in analyzing broad-band data of blazar sources constitutes a very 
promising method for investigating particle acceleration processes operating in relativistic plasma. 
In the particular case of the BL Lac object Mrk 421, it already allowed us to limit the number of 
possible scenarios to be considered in this context. Still, several different options remain, which 
could be investigated in detail using the SYNCHROTRON model with new broad-band, good-quality 
UV-to-X-ray data covering long time intervals and different activity states in this and other sim- 
ilar sources. In this paper we focused on the simultaneous UV to X-ray Mrk 421 data from Swift 
and Suzaku observations in December 2008, and we argued that gamma-ray data simultaneous to 
the UV- X-ray data would add further constrains to the various possible model/physical interpre- 
tations. We want to point out that there exists Fermi/MAGIC observations contemporaneous to 
the Swift/ Suzaku observations during December 2008. The analysis of this data set, as well as 
the physical interpretation of the overall SED in the context of the above presented model will be 
reported on a forthcoming publication. 

L.S. and M.O. are grateful for the support from Polish MNiSW through the grant N-N203- 
380336. GM acknowledges support via NASA grants NNX08AX77G and NNX08AZ89G. 
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Fig. 1. — The UV-to-X-ray spectral energy distribution of Mrk 421 fitted by the SYNCHROTRON 
model. The top panel shows the analyzed data fitted with (a) a power-law EED with fixed 7 m i n = 1 
(solid curve), (b) a power-law EED with 'y m i n set free (dashed curve), and (c) a modified ultra- 
relativistic Maxwellian EED (dotted curve). The bottom panels represent the residuals between the 
data and the models. The XIS data points are represented by filled circles in order to distinguish 
from the XRT and UVOT data. 
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Fig. 2. — The dependence of the electron power-law slope (cases "a" and "b" represented by filled 
and open circles, respectively) on the Galactic column density, iVn (top panel). The extinction 
parameter E(B-V) is scaled linearly with A^h (see text). The allowed range for ./Vh inferred from 
the historical HI observations is indicated by a red rectangle. The corresponding reduced x 2 values 
are plotted in the bottom panel. 



